z0 = -1:0.1:3;
r = exp(-z0.^2/2) .* sin(z0);
[x, y, z] = cylinder(r);
z = -1 + 4 * z; 
surf(x, y, z) 